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Introduction 

The  goal  of  this  study  is  to  identify  novel  imaging  marker  such  as  diffusion  kurtosis  that  may  serve 
as  a  prognostic  marker  in  the  evaluation  and  management  of  traumatic  brain  injury  patients.  This 
study  will  recruit  patients  patients  with  varying  severity  of  traumatic  brain  injury  (mild,  moderate, 
severe)  from  the  Shock  Trauma  center  with  an  initial  diagnosis  of  diffuse  axonal  injury. 
Neuropsychological  testing  and  magnetic  resonance  imaging  will  be  performed  in  the  acute,  sub¬ 
acute  and  chronic  stages  (through  6  months).  The  relationship  between  the  advanced  magnetic 
resonance  imaging  markers  and  the  clinical  condition  of  the  patient  will  be  evaluated  at  each  time 
point  to  determine  which  of  the  imaging  markers,  or  a  combination  of  imaging  markers  are  best 
representative  of  the  clinical  condition  of  the  patient.  Further,  the  markers  will  be  evaluated  for 
their  prognostic  ability  to  determine  the  clinical  course  of  the  TBI  patients. 

Overall  Progress 

Significant  delays  were  experienced  in  the  establishment  of  the  award  and  the  approval  of  the  IRB 
protocol  for  this  project.  However  we  are  now  in  an  active  state  of  recruitment  of  patients  which  is 
progressing  smoothly. 

Protocol  Modification 


The  IRB  approvals  from  University  of  Maryland  and  eventually  from  the  HRPO  took  significantly 
longer  than  expected.  Some  of  the  delays  were  related  to  coordination  between  our  site  and 
Georgetown  which  unfortunately  has  been  in  a  flux.  Since  Dr.  Mhyre  leaving  Georgetown 
sometime  in  December  there  have  been  two  other  personnel  changes  which  delayed  the  start 
date  for  recruitment.  Currently  we  have  Dr.  Fiandaca  as  the  contact  person  for  receiving  the  blood 
samples  at  Georgetown.  This  modification  was  made  as  recently  as  March  15*'^,  2013.  We 
believe  we  have  all  the  approvals  now  for  us  to  actively  recruit  patients 

Study  Tasks 

Screening  and  recruitment  for  this  study  has  begun. 

Task  1 

a.  Enroll  patients  into  the  study.  Screening  and  recruitment  for  this  study  has  begun.  We  have 
already  recruited  16  patients  into  the  study,  but  unfortunately  we  were  not  able  to  bank  the  blood 
samples  from  these  patients  because  the  bio-sample  recipient  site  was  not  ready. 
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b.  Obtain  Diffusion  Kurtosis  images  within  10  days  of  injury,  4  weeks  and  6  months 
foliowing  injury.  For  the  enrolled  patients,  each  has  completed  an  initial  MRI  (within  10  day).  For 
some  of  these  patients,  only  a  standard  of  care  MRI  was  able  to  be  obtained  for  various  clinical 
reasons,  but  for  the  majority  we  have  obtained  DKI  images.  Several  of  the  enrolled  patients  have 
returned  for  their  4  week  follow-up  visit  and  undergone  DKI  scans.  We  will  continue  to  follow 
enrolled  patients  to  obtain  follow-up  scans. 

c.  Obtain  data  on  normal  control  subjects.  3  healthy  volunteers  have  been  enrolled  as  normal 
control  subjects,  and  DKI  images  have  been  obtained. 

d.  Obtain  blood  samples  at  the  time  of  each  MRI.  Due  to  the  above  mentioned  delays  in 
coordinating  with  our  bio-sample  recipient  site,  we  have  not  obtained  blood  samples  from  the 
enrolled  patients.  However,  we  are  now  in  a  position  to  begin  collecting  blood  samples  from 
prospectively  enrolled  patients. 

e.  Pre-process  the  imaging  data.  Convert  ail  imaging  data  to  standard  atlas,  perform 
regional  analysis,  perform  histogram  analysis  for  whoie  brain  data,  quantify  DTI  and 
kurtosis  anaiysis  quantify  MRSI  data,  deveiop  hemorrhagic  iesion  burden  scale,  perform 
image  segmentation  to  separate  white  and  gray  matter  and  obtain  regional  volumes.  Pre¬ 
processing  has  begun.  As  imaging  data  begin  to  accumulate  we  will  begin  to  perform  various 
analyses. 

Processing  of  data  will  continue  on  an  ongoing  basis  for  both  data  quality  monitoring  purposes 
and  to  assess  the  general  trend  of  the  data  as  new  data  is  obtained.  Data  is  reviewed  every 
Thursday  afternoon  with  the  team  to  discuss  various  ways  to  analyze  data  and  potential  for 
manuscript  generation. 

A  technical  presentation  on  the  behavior  of  DKI  when  multiple  diffusion  weightings  and  multiple 
directions  are  used  was  presented  at  the  International  Society  for  Magnetic  Resonance  in 
Medicine.  We  investigated  whether  the  number  of  directions  were  more  important  compared  to 
the  number  of  diffusion  weightings  (b-factors)  in  order  to  sample  the  kurtosis  tensor.  Our  interest 
was  also  to  determine  the  shortest  acquisition  time  possible  without  losing  much  information 
regarding  the  tensor.  It  was  determined  that  the  maximum  b-value  of  2500  s/mm^  and  a  minimum 
of  two  b-values  provided  a  reasonable  representation  of  the  tensor  for  an  acquisition  of  time  of 
less  than  ten  minutes.  The  slides  from  this  presentation  are  provided  in  the  Appendix.  A 
manuscript  of  this  work  is  near  completion  for  submission  to  Magnetic  Resonance  in  Medicine 
Journal. 


Work  Durinp  Next  Reporting  Period 

We  will  now  actively  recruit  patients  for  the  study.  We  expect  to  recruit  at  the  rate  of  about  10-15 
patients  a  month  which  is  about  the  average  from  the  ongoing  parent  study.  The  parent  study  will 
stop  recruiting  any  new  patients  and  all  the  new  patients  will  flow  into  this  study.  By  the  next 
quarter  we  hope  to  have  about  40-50  patients  recruited  with  about  half  the  patients  finishing  their 
one  month  follow  up. 

Administrative  Comments:  We  hope  that  all  the  administrative  issues  are  behind  us  and  can 
actively  recruit  patients.  At  this  time  there  are  no  administrative  comments. 
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Introduction 


•  Diffusion  Kurtosis  Imaging  (DKI)*  measures  the  non- 
Gaussian  behavior  of  water  diffusion  and  has  gained 
much  interest  lately  as  a  tool  to  reveal  significantly  more 
tissue  microstructure  information,  over  and  above  which 
can  be  achieved  from  than  Diffusion  Tensor  Imaging 
(DTI) 

•  However  clinical  application  of  DKI  faces  a  major 
challenge,  as  it  involves  long  acquisition  times  (~20  min) 
due  to  more  complex  model  (21  model  parameters  vs.  6 
in  DTI) 

•  Little  is  known  about  how  different  DKI  imaging 
schemes  affect  the  variability  in  estimated  DKI 
parameters 


*  Jensen  JH,  et  al.  2005 
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Goal  of  this  study 

•  To  study  the  effect  of  imaging  schemes: 

°  Diffusion  weightings  (b-values) 

°  Diffusion  directions 
on  estimated  DKI  parameters 

•  To  find  an  optimal  DKI  imaging  schemes  within 
a  clinically  feasible  image  acquisition  time 

(<  lomin) 
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Diffusion  Weighted  Imaging  Schemes 


Full  Dataset 


4  DKI  dataset,  each  with:  b2  b^  b^  b^ 

5  b-values:  b  =  500, 1000, 1500,  2000,  2500  s/mm^ 
64  diffusion  directions 

Acquisition  time:  30  minx 4 


Diffusion  Weightings 
(b-value) 


Diffusion  Directions 


Optimal  imaging  schemes 


Diffusion  weighted  imaging  schemes  were  constructed 
by  taking  subsets  from  the  full  DKI  datasets. 

Experiment  1:  Find  the  optimal  b-value  subset 

Experiment  2:  Find  the  optimal  combinations  of 
numbers  of  b-values  and  diffusion  directions 
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1 

Methods 


•  Both  regions  of  interests  (ROI)  and  whole  brain  analysis 

•  For  ROI  analysis,  four  regions  are  selected: 


^  rm. 


Genu 


Internal  capsule  Thalamus 


{G  V  T  ^  I-  € 

X'  it  >  _n  ‘  i. 

•  m' 


¥ 


Basal  ganglia 


•  In  order  to  minimize  directional  bias  for  the  chosen  Ndir 
sets  to  the  tensor  direction  in  each  ROI,  eight  maximally 
different  subsets  of  15,  30  and  45  diffusion  directions  were 
chosen  from  the  fall  64  direction  set 

-Based  on  minimizing  electrostatic  energy  of  the  direction  vectors 
-  32  data  points  for  each  ROI  (4  repeats  x  8  subsets) 
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Effect  of  Diffusion  Weighting  Selections 


6.00  n 


0.00 
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NbvaL 
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Nbvai^  Nbvalg 


Estimation  Variability 


NbvaL 


NbvaL 


Nbval^  Nbvalg 


•  All  sets  with  of  2000  s/mm^  have  higher  bias  and  variability 

•  Optimal  b-value  subsets: 

Nbvalg:  h  =  1000,  2500  s/mm^ 

Nbvalg:  h  =  1000,  2000,  2500  s/mm^ 

Nbval4:  h  =  1000, 1500,  2000,  2500  s/mm^ 
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N  I  I  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1 


I  1  1  1  1  1  1  1  1  I  1 

1  1  1  1  1  1  1  1  1  1  I  I 


Effect  of  Numbers  of  b-values  and  Diffusion  Directions 

Genu  Thalamus 


o\ 

5 

QJ 

bJO 

cd 

(D 


4.50  - 
4.00  - 

3.50  - 
3.00  - 

2.50  - 
2.00  - 

1.50  - 
1.00  - 
0.50  - 
0.00  - 


*'  V* 

I 


14.00 

12.00 

10.00 

8.00 

6.00 

4.00 

2.00 

0.00 


Nbval2  “►  Nbvalj 


♦  Ndir, 


64 


♦  Ndir 


45 


♦  Ndir 


30 


♦  Ndir 


®(j)  V*  ♦♦  ♦  ♦ 


15 


O 


10 


20 


30 


40 


10 


20 


30 


40 


Imaging  acquisition  time  (min) 


Imaging  acquisition  time  (min) 


Opt7min:  Nbval  =  2(b  =  1000,  2500  s/mm^),  Ndir  =  30 
Optiomin:  Nbval  =  2  (b  =  1000,  2500  s/mm^),  Ndir  =  45 


•  How  do  they  perform  compared  to  already  proposed  DKI  imaging  schemes? 

Origiymin^:  Nbval  =  5  (b  =  500, 1000, 1500,  2000,  2500  s/mm^),  Ndir  =  30  (-17 
Orig7min^:  Nbval  =  2  (b  =  1000,  2000  s/mm^),  Ndir  =  30  (-7min) 

Gold  Standard:  Full  data  (Nbval  =  5,  Ndir  =  64,  average  of  4  repeats) 

^  Jensen  JH,  et  al.  2005;  ^Jensen  JH  &  Helpern  JA.  2010.  7 
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•  Optiomin  >  Optymin  «  Origiymin  >  Origymin 
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Estimation  Accuracy  in  Whole  Brain 
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•  The  Origymin  scheme  has  much  higher  bias  due  to  a  different  (not  shown) 

•  Accuracy:  Origiymin  >  Optiomin  >  Optymin 

•  Bias  for  Optiomin  may  be  systematic  (lower  spatial  variability) 

•  Negative  bias  in  MK  and  is  noise  related 
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How  Does  Image  Noise  Affect  DKI  Fitting 

•  Diffusion  parameters  have  to  be  physically  meaningful. 


Principle  of  diffusion  MRI: 

>^0  ^  Constraint  i:  D>0 


HS(b)/s,) 


Constraint  violations  for  DTI  parameters  are  Rare! 


D<0 


D>0 


>Sib,)>S(b, lift, <b, 
Constraint  2: 


HS{b)/S,) 


b 


From  multi-compartment  model  and  empirical  evidence*: 

>  \n(S(bySQ)  is  a  convex  function  of  b(K>  0)  in(^(Z))/5’o) 


=>  Constraint  3:  K  >0 


K>0 


Constraint  violations  for  DKI  parameters  are  very  common! 


t  K<0 


b 


*Tabesh  et  al.,  Magn  Reson  Med.  2011 
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Constraint  Violations  of  Different  Imaging  Schemes 

Voxels  with  K<0 

Method:  Count  voxels  which 
violated  different  constraints* 

-  Negative  diffusivity  (D  <  o) 

-  Negative  kurtosis  (K  <  o) 

-  Kurtosis  value  over  limit 
(K< 
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•  Voxels  with  i)<  0  is  rare  (  <  0.5%) 
(now  shown  here) 


•  Percent  voxels  violations: 
Optiomin  <  Origiymin  < 
Optymin 


WM 


CSF 


*Tabesh  et  al.,  Magn  Reson  Med.  2011 
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Conclusion 

•  Overall  DKI  estimation  benefit  more  from  increased 
number  of  diffusion  directions  than  number  of 
diffusion  weightings. 

°  More  diffusion  directions  reduces  estimation 
variability  of  DKI  parameters. 

°  More  diffusion  weightings  increases  estimation 
accuracy. 

•  The  optimal  DKI  imaging  schemes  are  with  two 
diffusion  weightings.  With  45  diffusion  directions, 
the  optimal  imaging  scheme  even  improves  upon 
the  conventional  imaging  with  17  min. 
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